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Abstract 



A new mechanism of the Verway transition in magnetite (Fe304), which has 
been argued to be a charge ordering transition so far, is proposed. Based on 
the mean field calculations for the three band model of spinless fermions ap- 
propriate for the d electrons of the Fe ions on the B sites, it is indicated that 
the phase transition should be the bond dimerization due to the cooperative 
effects of strong electronic correlation and electron-phonon interaction. The 
results show that the ferro-orbital ordered state is stabilized in a wide temper- 
ature range due to the strong on-site Coulomb interaction between different 
t2 g orbitals resulting in an effectively one-dimensional electronic state, which 
leads the system toward an insulating state by the Peierls lattice distortion 
with the period of two Fe(£>) ions, i.e., bond dimerization. Furthermore, it 
is found that the interplay between such lattice distortion in the Fe(B) ions 
and the lattice elastic energy of the Fe(B)-0 as well as the Fe(A)-0 bonds 
gives rise to a competition between two different three-dimensional patterns 
for the bond dimerization, and can stabilize a complicated one with a large 
unit cell size. The results are compared with the known experimental facts. 

PACS: 71.10.Fd,71.30.+h,71.10.Hf,75.50.Gg 
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I. INTRODUCTION 



The nature of the Verway transition in magnetite (FesC^) has not been clarified yet, in 
spite of intensive studies from very early days ||. Fe 3 4 forms the cubic spinel structure, as 
shown in Fig. [I], where one-third of the Fe ions occupy the A sites, tetrahedrally coordinated 
by 4 oxygen ions, while the remaining two-third are octahedrally surrounded by 6 oxygen 
ions, which are the B sites. One can see that Fe(5)204 layers and Fe(A) layers are stacked 
alternatively, while the z axis can be taken as any of the three cubic crystal axes. Although 
the primitive cell is a rhombohedral parallelepiped containing 2 formula units of Fe3C>4, the 
unit cell is conveniently chosen as the cubic one, as shown in Fig. [I], containing 8 formula 
units, to which we also refer in this paper. 

The Fe(A) ions are trivalent while the Fe(B) ions are mixed valent with formal average 
valence Fe(5) 2 ' 5+ . Below Tn = 858 K, the magnetic moments of the Fe ions are ferrimagnet- 
ically ordered, where the A sites and the B sites have opposite spin directions, with <i-orbital 
occupation represented as (t2 3 t) 3 ( e pT) 2 anc ^ (hgl) 3 '( e glY ! (^2 9 t)°' 5 ' respectively. Well below T N , 
an abrupt increase in the resistivity takes place at Ty — 120 K, which is now called after 
Verway, who has discovered it |2| and proposed that the phenomena is due to the ordering of 
equal number of Fe 2+ and Fe 3+ on the B sites, i.e., charge ordering (CO) among the "extra" 
t 2g i electrons |J. 

As for theoretical studies on this phenomena, the electronic properties of the spin po- 
larized t 2g ] electrons have been investigated frequently by the single band spinless fermion 
model on the B sites originally proposed by Cullen and Callen f|||. Such a picture of 
spinless fermion has been proved to be valid by band structure calculations , where the 
majority spin (J,) band of the Fe(B) ions is shown to be fully occupied and the Fermi energy 
crosses the minority spin ( j) band. In the model of Cullen and Callen the triple degeneracy 
of the ti g orbitals is neglected so that the band is half-filled, corresponding to half a charge 
per site as in the actual compound. This model is investigated within the Hartree approxi- 
mation where the ground state shows a phase transition from a metallic state to a CO state 
as the intersite Coulomb interaction V%j is increased [|J , providing the most naive picture for 
the CO state in the magnetite. 

However, no CO pattern so far considered has been successful in presenting a satisfac- 
tory interpretation for the known experiments, including the one originally proposed by 
Verway || and other more complicated patterns |]9|JT0[|. Besides, the expected short range 
fluctuation of the CO due to the frustration in Vij [11], arising from the fact that the network 



connecting the Fe(B) ions forms a coupled tetrahedra system, is not found either in neutron 
scattering [ll| or in resonant X-ray scattering measurements [ i~3|[ . Furthermore, the recent 
NMR |nj and X-ray anomalous scattering [ffH] experiments for T < Ty even cast doubts 



on the existence of CO. Todo et al. ||16|| , based on these suggestions and their finding of 
pressure- induced metallic ground state above 8 GPa, proposed that the low temperature 
phase below Ty may be a kind of "Mott insulator" where the B sites are forming dimers 
due to orbital ordering (00). 

This idea of the Mott insulating state resembles to that in low-dimensional quarter-filled 
organic conductors theoretically studied by Kino and the present authors Hl7j , where the 



fact that there exists a carrier per 2 sites is the same as that for the B sites in Fes04. In 
such organic conductors, the Mott insulating state can be understood from a view point 
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FIG. 1. Cubic spinel structure of Fe3C>4. The thick square shows the cubic unit cell with a 
lattice constant a in the xy plane while the z coordinates are indicated by the fraction of a. The 
unit cell in the xy plane for T < Ty is shown by the dotted line. 



that each charge localizes in every two sites, i.e., dimer. 



to the anisotropy in the transfer integrals from the lattice structure [18 



The dimerization is either due 
or due to the 



spontaneous formation of bond dimerization (BD) by the electron-phonon interaction where 
the one-dimensionality (ID) plays a crucial role [|19| . 

In this paper, we will show that this latter type of insulating state with BD can actually 
emerge in Fes04, and propose it to be the mechanism of the Verway transition in this 
compound. We will discuss in Sec. ||that the strong correlation among electrons stabilizes an 
00 state with effectively ID electronic structure, and additional electron-phonon interaction 
can give rise to the BD state. Moreover, discussions in Sec. |T| on the elastic lattice energy 
of the whole system will lead to the three-dimensional pattern of this BD expected below 
Ty. The relevance of our proposal to the experimental facts is discussed in Sec. |TV], and the 
conclusion is given in Sec. |V|. 



II. BOND DIMERIZATION ON B SITES 
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A. Three band model 



We start with the model of Mishra et al. [pOl , who have extended the spinless fermion 
model of Cullen and Callen j3J by including the triple degeneracy of the ti a orbitals appro- 
priate for the B sites in Fes04. The Hamiltonian is written as 

+ J2Vn inj , (1) 

(ij) 

where cj (c ifJ ) and n iM denote the creation (annihilation) and the number operator of the 



electron at zth site of orbital fi, respectively, where the orbital indices take xy, yz or zx [FIT 
rii is the number operator for each site, i.e., rii = n ixy + n iyz + rii ZX , and (ij) denotes the 
nearest-neighbor site pair along the coupled Fe(B) ion tetrahedra network. U and V are 
the on-site Coulomb energy between different orbitals (note that the same orbital cannot 
be doubly occupied) and the nearest-neighbor Coulomb energy, respectively. By using the 
transfer integrals for (ij) pairs calculated in Ref. [7], tdda = —0.41 eV, tddn = 0.054 eV and 
tdds = 0.122 eV, the three kinds of transfer integrals, t^, with different configurations of 
orbital occupations, as shown in Fig. |2], can be estimated as t\ = 3tdda/^ + tdds/^ — 
-0.278 eV, t 2 = W2 + W2 = 0.085 eV and t 3 = t dda /2 - t dd5 /2 = -0.035 eV, where the 
notations of t lr ^ 3 are given in Fig. [2|. 

Mishra et al. treated the Coulomb interaction terms, U and V, by means of Hartree 
mean field approximation, 

n^riju -> (riifj) n jv + n ifl (n ju ) - (n if ,) (n ju ) , (2) 

and determined the ground state of the system. Their results for the case with the fixed 
value of U = 4.0 eV, relevant for magnetite 0, showed that the CO state is stabilized 
in the ground state for V > V c = 0.38 eV |2^| while the metallic state is for V < V c , 
qualitatively the same as in the one-band model M, and they concluded that the former 
state is relevant for the magnetite. 

Although it is not emphasized in their paper, every site is almost fully occupied by the 
same orbital state, namely, ferro-00 is realized in both of these states, due to the large 
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FIG. 2. A schematic representation of the transfer integrals between neighboring B site pairs in 
the xy plane with different orbital occupations. The simultaneous cyclic substitution of the orbital 
indices and the axes of coordinates, x — ► y, y — > z, z — ► x or x — ► z, y — ► x, z — ► y, also provides the 
three transfer integrals, while for the other configurations t^=0. 
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value of U. The three d xy -, d yx - and d yz -00 states are energetically degenerate, and the 
electronic structure becomes ID once one of the 00 state is realized, which is particular to 
the spinel structure. For example the d xy -00 state has ID arrays of sites connected by the 
largest transfer integral t\, which are along the [110] direction and the [110] direction for the 
xy planes with z coordinates z = (4n + l)/8 and z = (4n + 3)/8, respectively (see Fig. [I]). 

In the metallic state with the presence of this d xy -00, the expectation values of the 
charge density for all sites are (rii xy ) ~ 0.5 and (n iyz ) = (n izx ) ~ 0. As for the CO state, i.e., 
the coexistent state of OO and CO, the planes with ID arrays of (n ixy ) ~ 0.5 + 5 and those 
with (n ixy ) ~ 0.5 — 5, corresponding to Fe 2+ and Fe 3+ , respectively, are stacked alternatively 
along the z direction, where the CO pattern is the one proposed by Verway. Here 5 is the 
amount of charge disproportionation, which rapidly increases as the value of V is increased 
from V c , and reaches around 0.9 for V = 1.0 eV. The other two d yz and d zx orbitals are 
almost empty also in the CO state. 

In contrast to their conclusion that this CO state is realized in magnetite, the recent 
experimental proposals of the absence of CO mentioned in Sec. | suggest that the former 
OO metallic state should be relevant to the actual system, and that the origin of insulator 
is other than CO. The destabilization of the CO state may be due to the screening of 
the long range Coulomb interaction and/or the effect of frustration among mentioned 



above [11 1. Actually, the melting of CO due to frustration among has been demonstrated 



theoretically in ID systems [24], though we will not discuss such possibility in magnetite 
further in this paper. Here we will concentrate in the following on how the OO metallic 
state shows instability toward an insulating state and its possible consequence on the lattice 
structure in this insulating phase. 



B. Peierls instability in orbital ordered state 

In the following we consider the case of the d xy -00 metallic state. The properties of this 
ferro-00 state can be extracted by an effective noninteracting ID spinless fermion system, 
ifiD = J2iti{cjci + i + h.c), with a half-filled band, / being the site index along the ID 
directions [110] or [110] . This model should be valid for U t\ t 2 ,t 3 in eq. (fj) with 
Vij = 0. It is well known that such a half-filled ID band has the Peierls instability with the 
wavelength of two times the inter-atomic distance ||19|| , i.e., instability toward the BD state 
with alternating atomic displacements along the ID direction. Actually the system can be 
mapped onto the 5=1/2 XY chain via Jordan- Wigner transformation, so that our metallic 
state here correspond to the spin liquid state, which shows instability toward the spin- Peierls 
state, that is, the BD state in spin system. Thus, once the electron-phonon interaction is 
present, the OO state in magnetite will undergo a BD transition to gain kinetic energy by 
making a gap at the band center, resulting in an insulating state. We propose this scenario 
to be the origin of the insulating ground state of magnetite, Fe 3 4 . Such a state will be 
realized if it is not destroyed by the inter-chain interactions t% and t 3 , which is investigated 
in the next subsection. 

We should note that this Peierls instability is different from the usual one in the ID 
electron system of weak coupling, since the assumption of the spinless fermion here is due to 
the ferrimagnetic spin ordering realized in the limit of infinite value of on-site intra-orbital 
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FIG. 3. Mean field phase diagram for the three-band Peierls-Hubbard model for the B sites, 
on the plane of temperature, T, and the inverse of the lattice elastic coupling constant, 1/K. OO 
and BD represent orbital-ordering and bond dimerization, respectively. 



Coulomb interaction. Thus one may say that the BD state here is rather analogous to the 
Mott insulator in quarter-filled compounds noted in Sec. H ||i~7Hi9|l, where the origin is also 
the interplay between the strong on-site Coulomb interaction and the dimerization. 



C. Peierls-Hubbard model 

To investigate the stability of the BD state in the actual three-dimensional system, we 
include the lattice degree of freedom by adding the Peierls-type coupling to Eq. (|l|), i.e., 
treat the Peierls-Hubbard model for magnetite. It is expressed as 

#ph = EEC( 1 + u ij) c tcju + E Y, unipniv 

<ij> V-v i 

+ E ~^4> ( 3 ) 

i 

where Uij and K are the lattice distortion between nearest-neighbor site pair (ij) and the 
coupling constant for the elastic energy, respectively. The intersite Coulomb interaction 
term is neglected here since it is not relevant in our discussion as mentioned in Sec. flA . 



The on-site Coulomb interaction term, U, is treated within the mean field approximation, 
as in eq. (§), and self-consistent solutions are obtained. We restrict ourselves to the case 
of = (— 1) u associated with the three kinds of transfer integrals along the [110] and 
[110] directions and u^j = for the other directions. This provides the BD state, u being 
the degree of BD, and the value of u is determined for each choice of parameters so as to 
minimize the energy. We note that the calculated energy does not depend on different BD 
patterns since the effect of the lattice distortion on the lattice elastic energy and on the 
transfer integrals is restricted in each chains. 

The results show that the BD state can actually be stabilized when 1/K exceeds a 
critical value. In Fig. [| the obtained phase diagram within the finite temperature mean 
field approximation is shown for the case of fixed U = 4.0 eV, the same as in the calculation 
of Mishra et al. [20[ mentioned in Sec. |11 A|. The ground state is OO metallic with u = 
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FIG. 4. Temperature dependences of (a) the degree of BD, u, and (b) the band gap for several 
values of 1/K. 



below 1/K < 1.9 eV -1 , whereas the 00 state with u ^ 0, i.e., the BD state, is stabilized 
for 1/K > 1.9 eV _1 , where the system is insulating. Below Too — 6000 K, which is the 
temperature range shown in Fig. ^ the ferro-00 is present, and the para-orbital state is 
stabilized only above T G o- 

The optimized value of u and the band gap in the BD state are plotted as a function of T 
in Figs. [|(a) and^(b), respectively. It can be seen that the lattice dimerization and the band 
gap decrease as the temperature is increased and vanish continuously at T = T c , showing a 
second order insulator-to-metal phase transition. On the other hand, the charge density for 
each orbital state does not change noticeably from (n ixy ) ~ 0.5 and (n iyz ) = (n izx ) ~ 0, for 
all sites, in the temperature range we are concerned. In other words, the OO is not affected 
through the metal-insulator transition. 



III. BOND DIMERIZATION PATTERN 

The three-dimensional BD pattern in the actual compound cannot be determined by the 
preceding calculation alone, since the calculated energy does not depend on the inter-chain 
configuration, as mentioned above. To discuss theoretically the stability of different BD 
states, the lattice elastic energy should be the most important factor since the BD produces 
large lattice distortion so that it moderately affects the lattice energy. The influence of the 
BD on the the inter-chain transfer integrals would be small so the difference in the kinetic 
energy between different BD states is neglected in the following discussion. 

The lattice elastic energy in Fes04, £i at , can be estimated by the sum of y(A«) 2 for all 
the nearest neighbor Fe-0 bonds, where Au is the deviation of the Fe-0 distance from that 
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in the equilibrium position above Ty and K is the elastic constant. The value of K should 
take common values K A for the Fe{A)-0 bonds and K^ for the Fe(S)-0 bonds, since all the 
Fe(^4)-0 bonds as well as the Fe(B)-0 bonds are respectively crystallographically equivalent 
above Ty. Thus £[ at can be expressed as £ A iat + ^Biat, for the Fe(A)-0 and the Fe(B)-0 
bonds, respectively, where 



^Alat 
^Blat 



E^(An A _ ) 2 , 
E^(A MB -o) 2 , 



(4) 
(5) 



where Awa-o and ^ u b-o are Am for the Fe(v4)-0 and Fe(B)-0 bonds, respectively. 

Then, as will be explained in detail later, there are two candidates for the BD states 
costing much less elastic energy than the others, as schematically shown in Fig. |5|(a) and 
(b). In the former the dimerization pattern between adjacent ID chains in the xy plane is in- 
phase, thus we call it the in-phase BD state, where the unit cell is unchanged from the cubic 
unit cell above Ty. On the other hand, in the latter pattern the BD is anti-phase, which is 
called anti-phase BD state in the following. Here, the unit cell becomes large as a/2 x y/2 x 2 
of the cubic one, as shown in the Figure. We will see below that the competition between 
two BD states arises depending on the relative value of K& and Kb, and our proposal is 
that the anti-phase BD pattern is realized in the actual compound. 

To see this competition, let us calculate a semi-phenomenological Landau-type free en- 
ergy by taking into account of both in-phase and anti-phase BD states. We consider the 
amount of lattice distortion for the BD in the Fe(B) ions to be uniform along each chains, 
denoted by u% and u 2 , for alternative chains in the xy planes. Then the in-phase BD is 
characterized by u\ = u 2 while the anti-phase BD is by u% = —u 2 , as shown in Fig. |^(a) 
and (b), respectively. The motions of oxygens are approximated to be perpendicular to the 
chains in the xy planes, parametrized by uq\ and uo 2 , for the chains with BD of u\ and u 2 , 
respectively. Within these approximations, the Fe(A) ions becomes all crystallograhpically 
equivalent and we allow their motions in any direction, thus represented by the x, y, and z 
components of the displacement vector, u Ax , u Ay , and u Az , respectively (see Fig. H). 

Then the free energy per formula unit of Fe3C>4, T , can be computed, where there are 4 
Fe(A)-0 bonds for ^Aiat and 12 Fe(5)-0 bonds for £siat, as 



J~ — ^BD + ^Alat + ^Blat; 

J'bd = az(T - T c )(u\ + u\), 



(6) 
(7) 



F Ka 

C-Alat — ~TT 




-^+U A y~ Uqi J + + It A; 



.V2 



UAy + U 02 



UK, 



—7=-UAx + U Q1 \ + U A y + I - - U A 
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(a) in-phase BD state 
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(b) anti-phase BD state 

z=l/16 z=9/16 



z=3/16 



z=5/16 





z=ll/16 



z=13/16 



z=7/16 z=15/16 




FIG. 5. Schematic representation of (a) the in-phase BD state and (b) the anti-phase BD 
state both coexisting with the d xy -00. Only the Fe(l?)204 layers are shown, where the thick 
Fe(-B)-Fe(-B) bonds represent the dimers. The thick squares show the unit cell in the xy plane and 
the z coordinates indicate the fraction of unit cell size along the z direction. The unit cell remains 
unchanged from the cubic one above Ty shown in Fig. |l] for the in-phase BD pattern, while it 
becomes \/2 X y/2 x 2 of that for the anti-phase BD pattern. 
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(a) in-phase BD (b) anti-phase BD 




FIG. 6. Schematic representation of the ionic displacements in (a) the in-phase BD pattern 
and (b) the anti-phase BD pattern. Only two Fe(-B)20"3 and Fe(A) layers each are shown, and 
the directions of displacement are represented by arrows. The motions of the A sites in (a) is not 
shown, which are along the z direction. 
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where the lattice constant of the cubic unit cell is set to 4 so that the length of ¥e(A)-0 
and Fe(B)-0 bond without any distortion are 3/4 and 1, respectively. JF BD describes the 
instability toward the BD state along each chain as the temperature decreases, discussed in 
the previous Section. 

In Fig. 0, T as a function of (ui,u 2 ) is plotted for the optimized positions of the Fe(A) 
and O ions, obtained by minimizing it numerically by varying other variables, ua x , u^y, 
uaz, Uoi, and uo 2 - It is plotted for two sets of elastic constants (a) = 0.1, Kb = 0.3 
and (b) K& = 0.5, K B = 0.3, at several temperatures. Above T c , the shape of the parabolic 
curvature gets deeper as T c is approached, almost symmetrically along u\ = u 2 and U\ = —u 2 
suggesting that the fluctuations of both in-phase and anti-phase BD states develop. Below 
T c , for the parameters in Fig. |7|(a), the minimum of the free energy appears along u\ = u 2 , 
that is, the in-phase BD state is stabilized, while in Fig. ^(b) the mimimum is along u\ = —u 2 
which means that the anti-phase BD state is realized. Note that the phase transition is 
a second order one. The stability of these BD patterns can be understood as follows. 



In the presence of the in-phase BD pattern, all the Fe(-B)404 cubes [2S], connecting the 
Fe(i?)204 xy planes show ionic displacements as shown in Fig. ||(a), where all the interlayer 
Fe(-B)-0 pairs move respectively in the same direction, providing small Amb-o's. Thus the 
lattice elastic energy of the Fe(S)-0 bonds, ^Biat, should be the lowest for the in-phase BD 
states with such configuration, compared to other BD patterns containing cubes as shown in 
Figs. P(b) and |8|(c), with both Fe(-B) pairs forming dimers or both being inter-dimer Fe(-B) 
pairs, respectively. For example, the anti-phase BD shows all of these three kinds of cubes, 
realized with the ratio of 2 : 1 : 1, in order, thus costing rather higher £siat than in-phase 
BD. On the other hand, the Fe(v4)-0 tetrahedra would be quite deformed in the presence 
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(a) K A = 0.\,K B = 0.3 (b) K A = 0.5, K B = 0.3 
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FIG. 7. Calculated free energy, T, for (a) K A = 0.1, Kb = 0.3 and (b) K A = 0.5, K B = 0.3 
at several temperatures, as a function of the degree of dimerization for adjacent chains, (111,112). 
Note that the contour plot in the base is guide for eyes. 



(a) (b) 



Fe(B) 



O 




FIG. 8. Schematic view of the lattice displacements in the Fe(-B)-0 cubes connecting adjacent 
Fe(-B)204 layers when the BD state is stabilized. In the in-phase BD pattern only the displacement 
pattern (a) is realized, while in the anti-phase BD pattern (a), (b) and (c) are realized with the 
ratio of 2 : 1 : 1. 
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(a) in-phase BD 



O 



(b) anti-phase BD 



FIG. 9. Schematic view of the displacements in the Fe(A)-0 tetrahedra for (a) the in-phase 
BD pattern and (b) the anti-phase BD pattern. 



of the in-phase BD, which is shown in Fig. |9|(a), providing large lattice elastic energy of 
the Fe(A)-0 bonds, ^Aiat, while the anti-phase BD pattern make the ions surrounding 
the Fe(A) ions showing displacements as shown in Fig. [|(b), where the cost in <?Aiat will be 
rather smaller. This is because the motion of the oxygens are in the same direction for the 
latter pattern so that the Fe(A) ions can adjust their position to lower Aua-o, as can be 
seen in Fig. |](b). To summarize, one can say that <?Biat favors the in-phase BD state while 
^Aiat favors the anti-phase one. 

These naive discussions are consistent with the above calculation of the Landau-type free 
energy. There, we have observed that a competition arises between the in-phase and anti- 
phase BD states, stabilized for Ka ^ Kb and for > Kb, respectively. We propose that 
the latter is the case in magnetite, since in general, large deformation is hardly realized in 
the tetrahedra of oxygen surrounding the A sites of the spinel structure, which may be due 
to its tight packing compared to the B sites [^6|. The importance of such ionic displacement 



of the O ions surrounding the A sites in determining the electronic properties in the ground 
state has recently also been pointed out in another spinel compound AIV2O4 | 27fl . This 



anti-phase BD state provides the correct unit cell determined experimentally as will be 
discussed later. 



IV. COMPARISON WITH EXPERIMENTS 

Based on the discussions above, our proposal for the physical properties in magnetite 
Fe3C>4 can be summarized in Fig. [Hj which has the following new features: the existence 
of 00 in the B sites at all temperature range below room temperature, the BD fluctuation 
(both the in-phase and the anti-phase BD's) above the Verway transition temperarure, Ty, 
and the anti-phase BD stabilized in the insulating phase as a consequence of compensation 
with the lattice elastic energy of the whole system. These are compared with the known 
experimental facts in the following. 

Below the Verway transition temperature, the existence of ferro-00 is supported with 



dielectric measurements showing large anisotropy p9[ . Above Ty, in contrast, crystal struc- 
ture analysises show no evidence of the lowered symmetry from the cubic one, apparently 
contradicting with our prediction of the ferro-00 even in the metallic phase. However, 
the symmetry in the electronic structure has been pointed out to be lowered than the cu- 



bic one []30| based on magneto-crystalline anisotropy measurements |3T| and in the recent 
resonant X-ray scattering at room temperature [JT3], consistent with the 00 state. 

As for the BD, there is a strong support from neutron scattering measurement by 



Shapiro et al. 31^] , which is the observation of ID correlation above Ty, not explained 
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FIG. 10. Schematic phase diagram for magnetite Fe3C>4. The 00 and BD denotes the or- 
bital-ordered and bond dimerized states, respectively. The vertical axis is the temperature, T, 
while the horizontal axis can be interpreted either as the pressure, p, the "rigidness" of the lattice 
represented by the lattice constant, K, or the effective dimensionality. The bold line is the Verway 
transition temperature, Ty, varying as a function of such parameters. 



by the CO but consistent with the present BD picture. This ID nature of the BD transition 
can naturally provide an explanation for the critical fluctuation above Ty observed in large 
temperature range up to room temperature, such as in the pseudogap structure of the op- 
tical conductivity |£J and in neutron scattering measurements where it has actually been 
realized that the atomic displacement is playing an important role fll2|. Moreover, our sce- 



nario that the Verway transition is a phase transition from 00 metal to the coexistent state 
of 00 and BD explains the entropy change of AS* = _Rlog2 per mole of Fes04 estimated 
from specific heat measurements [34], arising from the degree of freedom for dimerization. 
This excludes the order-disorder-type CO picture of Verway H predicting AS = 2R log 2 as 
well as simultaneous 00 and BD transition with AS = 2i?log3 from the orbital degree of 
freedom. 

The crystal structure in the presence of anti-phase BD provides the correct unit cell size 
below Ty in the actual compound, namely, \/2 x y/2 x 2 of the cubic unit cell above Ty 
shown in Fig. [I] ||28|| , which is also supported by the double period modulation along the z 
axis observed in electron diffraction and in neutron scattering [ |36[| . Furthermore, the 
diffusive spot above Ty of the corresponding wave vector (0,0, 1/2), observed in a neutron 
scattering measurement showing a divergent behavior toward Ty jl2| can be assigned to the 
anti-phase BD fluctuation developping above Ty seen in the calculation of Sec. |HI . On the 
other hand, the analysis of Yamada et al. [[37] of their neutron scattering data above Ty 
apparently seems to correspond to the in-phase one, where the ionic motions of the Fe(B)-0 
cubes are proposed to be as in Fig. ||(a). This might be the observation of the fluctuation 
of the in-phase BD mode, which is also seen in our calculation though a long range order of 
this mode is not achieved. We have seen that such competition arises due to the existence of 
two kinds of elastic energies, ^Aiat and ^Biat, which may be consistent with the experimental 
findings by Shapiro et al. that "two types of interacting dynamical variables are taking part 
in the neutron scattering" [[32 



The observed pressure induced metallic behavior [16| can be naturally understood since 
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the pressure should effectively increase the value of K in Eq. (|]) (or Ka and Kb in Eqs. (£|) 
and (H) ) making the lattice more "rigid", and/or increase the three-dimensionality of the 
system, both expected to destabilize the BD state, which is represented in Fig. |lUj. This is 
analogous to the case of spin-Peierls transition, the mapped state of the BD state as men- 
tioned in Sec. |11B| , known to be destroyed by enhancement in the electron-lattice coupling 
constant and/or the three-dimensionality f38fl . We point out the possibility that the three- 
dimensional BD pattern may be changed under different environments such as pressure and 
temperature, since there are two competing phases, the in-phase and the anti-phase BD 
states, as seen in Sec. [TTJ. Search for such transition between these phases under uniaxial 
pressure is also interesting, since these states are highly anisotropic. Another possibility 
is that incommensurate phases may be stabilized as a consequence of such competition, as 
observed in dielectrics [N^CIDJ^MCU as well as in CO systems such as perovskite Ni 
oxides |4(J and NaV 2 5 fl41| , due to competing different states. 

One discrepancy between our theory and the experimental facts is the order of the phase 
transition: experimentally it is first order while it is second order in our calculations in 
Sees. || and [[II]. This discrepancy may be due to the simplified approximation we made for 
the strucural change as discussed in Sec. Il|. In the actual compound, the lattice distortion 
will be more complicated than in our calculation, represented as in Fig. |6|, e.g., the lattice 
distortion of the B sites along the chain should take the period of four sites which is contained 
in the unit cell, as (u a , —Ub,u c , —Ud,u a , —Ub, ...), rather than (u, —u,u, —u, ...) as in our 
calculation. This should make the crystal symmetry of the BD insulating phase lower than 
that in our calculation, which may lead to a first order phase transition j|2] . 

Finally, it is noted that the recent findings of possible charge disproportionation in the 
B sites by X-ray anomalous scattering ||3| does not contradict with our proposal here. The 
anti-phase BD state gives rise to crystallographically independent Fe(-B) ions, which should 
make the charge density on the B sites different where the amount of charge disproportiona- 
tion will be not so large but can be detectable in a sensitive probe such as X-ray anomalous 
scattering measurements. 



V. CONCLUSION 

In conclusion, we have theoretically proposed a new model for the mechanism of the 
Verway transition in magnetite FesO^ It is not a charge ordering transition as has been 
believed for more than a half century, but a bond dimerization induced by the Peierls 
instability in a one-dimensional state as a consequence of the ferro-orbital ordering due to 
strong electronic correlation. The actual bond dimerization pattern is predicted to be the 
anti-phase one, based on discussions on the lattice elastic energy of the system. This model 
seems to be able to provide the explanation for the long lasting mystery of the Verway 
transition in this compound. 
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